#include <stdio.h>
#include "itpack2c.h"

int main()
{
  int IA[5] = { 1,4,6,8,9};
  int JA[8] = {1,2,3,2,4,3,4,4};
  int IPARM[12], IWKSP[12];
  float A[8]  = {4.0, -1.0, -1.0, 4.0,
	       -1.0, 4.0, -1.0, 4.0};
  float RHS[4] = {6.0, 0.0, 0.0, 6.0};
  float U[4], WKSP[24], RPARM[12], zero=0.0;
  int N=4, NW=24, ITMAX=4, LEVEL=1, IDGTS=2, IER;
  
  dfault_(IPARM, RPARM);
  IPARM[0] = ITMAX;
  IPARM[1] = LEVEL;
  IPARM[11] = IDGTS;
  vfill_(&N, U, &zero);
  jcg_(&N, IA, JA, A, RHS, U, IWKSP, &NW, WKSP, IPARM, RPARM, &IER);
  
  return 0;
}
